Diffractive A+ Productions in Polarized pp 
Reactions and Polarized Gluon Distribution 



N. I. Kochelev [] 

Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 
141 980 Dubna, Moscow region, Russia 

and 

T. Morii f\ and S. Oyama | 

Faculty of Human Development and 
Graduate School of Science and Technology, 
Kobe University, Nada, Kobe 657-8501, Japan 



Abstract 

To test the model of the polarized gluon distribution AG(x, Q 2 ) in the proton, we propose a 
new process, diffractive A+ productions in polarized pp reactions, which will be observed in the 
forthcoming RHIC and also the proposed HERA-N experiments. The spin correlation between 
the target proton and the A+ produced in the target fragmentation region largely depends on 
AG(x, Q 2 ) and thus, the process is quite promising for testing the models of AG(x, Q 2 ). 
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In these years, spin structure of nucleons has been one of the most challenging topics in nuclear 
and particle physics]]]]. Since the surprising observation of the polarized structure function of the 
proton, gi(x), by the EMC collaboration in 1988, much progress has been attained theoretically 
and experimentally in the study of the spin structure of nucleons. A large amount of data 
on polarized structure functions of proton, neutron and deuteron were accumulated by many 
experimental groups such as SMC@ at CERN and E142@, E143||, E154§, E155[g at SLAC 
and HERMES [0] at DESY. The progress in the data precision is also remarkable. On the other 
hand, the next-to- leading order QCD calculations of polarized splitting functions || stimulated 
the theoretical activities for studying the polarized parton distributions in the nucleon. Several 
parameterization models of polarized parton distribution functions which fit well to the data have 
been proposed[|, [HJ [11], |12| and analyses have been developed with increasing new data|l3|, |i~4 



These experimental and theoretical developments brought about a deep understanding on the 
behavior of the polarized parton distribution in the proton. However, a knowledge of the polarized 
gluon distribution which is expected to play an important role in the so-called nucleon spin puzzle 
is still poor, though many processes have been proposed so far to extract information of it. 

In this work, in order to extract information on the polarized gluon distribution, AG(x, Q 2 ), 
we propose a different process, i.e. diffractive A+ productions in polarized pp collisions at high 
energies. Diffractive process which can be described by the Pomeron exchange is also an inter- 
esting current topic There have been many discussions on Pomeron interactions in these 



years, largely motivated by recent HERA experiments. To study this process is now quite timely 
because RHIC will start soon and one of the main purposes of RHIC experiment s ||16|| is to extract 
the polarized gluon distribution in the nucleon from various reactions sensitive to AG(x, Q 2 ) in 
the proton. We expect that RHIC will also observe our proposed reactions, the diffractive A+ 



production. The same process will be observed in the proposed HERA-N experiments! 17], too. 
The process which we are considering here is 

p(pO + Pfa) — > p(p[) + A+(p Ac ) + A, (1) 

where particles with arrows indicate that they are longitudinally polarized and parameters in 
parentheses denote the momenta of respective particles. The lowest order Feynman diagram for 
this process is shown in Fig.l. 

Let us start by describing why we are interested in A+ productions. It is well-known that 
the A+ is composed of a heavy quark c and antisymmetrically combined light u and d quarks, 



and thus, the spin of A+ is basically originated from the c quark. In addition, a c quark is 
produced from protons just through gluon fusion in the lowest order as shown in Fig.l and 
hence, gluon polarization affects the spin of A+ via c quarks produced from gluons. In other 
word, measurement of the spin of A+ gives us an information on the polarized gluon in the 
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The cross section of this process can be calculated based on the parton model by using the 
Pomeron model describing a diffractive mechanism and the model of polarized gluons, AG(x, Q 2 ), 
and polarized fragmentation functions, AD A + / c (£), of an outgoing c quark decaying into a polar- 
ized A+ with momentum fraction £. One of the conventional models of Pomeron interactions is 
the one proposed by Donnachie and Landshoff (DL) |19| , in which the Pomeron is considered to 
behave like a C = +1 isoscalar photon. Here we take this model as our first analysis. As a typical 
model of the polarized gluon distributions, we take GS96|| and GRSV96[Kj parameterizations 
among many models, since both of them reproduce well the experimental data on the polar- 
ized structure function of nucleons and thus, seem plausible. As for the polarized fragmentation 
function AD A +/ C (£), unfortunately it is not known at present because of lack of experimental 
data, though we have some knowledge of the unpolarized function -Da+/c(O0- However, since 
a c quark is heavy and hence, it is not expected to change much its spin alignment during the 
fragmentation process, it might not be unreasonable to use D A + / (£) for AD A + / c (£). Here we use 
the model of Peterson et al.[^(J as a substitute for AD A + / c (£). 



To calculate the cross section of this process based on the parton model, it is convenient to 
use the scaling variables, x = ~ A ! , y = E2 ^- and z = P2 ' Pt ^ where A = m — r>\ is the momentum 
transfer of the incoming unpolarized proton. The subprocess of this scattering is 

p(pi) + g(k) — ► p(p[) + c(p c ) + c(p g ), (2) 

where g and c(c) represent the gluon and c-quark(c-quark), respectively. Momenta of individual 
particles are given in parentheses and are related to the ones in the physical process of eq.(|]) 
as k = £p 2 , p c = ^"S where £ and £' are momentum fractions of the gluon to the target proton 
and the A+ to the c-quark, respectively. For this subprocess, we can also define the scaling 
variables, x p = = |, y p = |^ = y an( j Zp = |^ = ^. Here we follow the Schuler's 
way|ZT| to calculate the 3-body phase space for the subprocess of eq.(0). Then, we define the 
particle momenta in the final state cc quark(or Pomeron-gluon) C.M.S. in the subsystem (g), 
with p c + p 5 = A + k = 0, where the gluon momentum k and the target proton momentum p 2 



point to the positive z direction. An angle <fi between the proton plane (k x p\) and the c-quark 
plane (k x p c ) is defined by cos</> = ^^'}i x ^ c ) ■ In addition to s = (pi + p 2 ) 2 , it is convenient 
to define the Lorentz invariant kinematical variables, s = (k + A) 2 , t\ = (jpi — p[) 2 = A 2 and 
t 2 = (k — Pc) 2 , which for s ^> m 2 , can be expressed in terms of x, y, and z as s = (£ — x)ys, 
i\ = —xys and t 2 = —^yzs + m 2 , respectively. 

By using these variables, we can calculate the differential cross section for the process of 
eq.fllD as follows; 

dAa f x ™°* , f l d£ f 1 d£' f 27T ,. . „ ti . ^ dAa . ^ . . 

dx -T . -W dcf>AG(£,Q 2 ) - J J AD At/c (?), (3) 



dydz Jx min Juin £ -^Cin £' -'o dxdydz* 
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where the kinematical limits of the integral variables are given as x max = ^, x min = , 
^min = aj + ^ and £^ in = ^ (l - y/l - x max is determined from — 1 < ti which we 

take for ensuring the diffractive condition for the incident proton with momentum pi, x m j n from 
sin 7 > where 7 is the angle between the momentum of the incoming unpolarized proton and 
z axis, £ m i n from the condition s > 4m 2 and £' min from the requirement — 1 < cos 6 < 1 for the 
scattering angle 9 of the final c-quark in the subprocess (0), respectively. In eq.(§), the polarized 
subprocess cross section is given as 

dAa da ++ da^ da da-+ 

dx p dy p dz p d(f) dx p dy p dz p d(j) dx p dy p dz p d(p dx p dy p dz p d(p dx p dy p dz p d(l)'' 

- Vp \\M\ 2 -\M\ 2 +\M\ 2 __-\M l2 ^ 
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Vp (A\M Y \ 2 + A|7W 2 | 2 - 2Re{A(M l M* 2 )}), (4) 
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with 



-iM = -iM 1 -{-iM 2 ), (5) 
A\M li2 \ 2 = \M 1>2 \ 2 ++ - \Mw\l_ + \M 1>2 \ 2 __ - 1-^1,21-+, (6) 
A(M!M* 2 ) = (M 1 M* 2 ) ++ -(M 1 M* 2 ) + - + (M 1 M* 2 )—-(M 1 M* 2 )^+, (7) 

w ^ ere dx dy + dz dip an d 1-^1+-' f° r ex a m ple, denote that the helicity of the gluon and the c quark is 
positive and negative, respectively. A4\ is the amplitude corresponding to the subprocess shown 
in the Feynman diagram of Fig.l and explicitly written by 

[k - Pc) 2 - m 2 c 

o(i)-l 



/ o, \ a(t)—x 

xu(p[)3(3oF(t)f((k - p s ) 2 - m 2 ) 7 ^( Pl ) - , (8) 

\Sq/ 



where (3 and (3q are the charm quark-Pomeron and light quark-Pomeron couplings, respectively. 
s i = (p'i +Pc) 2 and So is a scaling constant fixed as Sq = lGeV 2 . M. 2 has similar expression which 
is originated from an interchange of c and c in the final state in the subprocess. Note that A4\ 
and are added with negative sign because of positive charge conjugation of the Pomeron, 
as shown in eq.(5). The unpolarized cross section can be calculated similarly by replacing the 
polarized functions in the integrand of eq.@, AG(^,Q 2 ), dx d ^ u dz and AD A + / c (£')> ^y the 
unpolarized ones, G(£,Q 2 ), dXpd ^ dZpd ^ and D A + /c (g), respectively 
The 2-spin asymmetry of A^" is calculated by 

_ da ++ — da + _ + da — da^ + dAcr da 

LL da ++ + da + - + da + da |_ dydz dydz 

By using these formulas, we have calculated the polarized and unpolarized cross sections and 
the 2-spin asymmetry for the diffractive A+ production in pp reactions at \fs = 50 GeV and 
y/s = 500 GeV, foreseeing the forthcoming RHIC experiments. 

Concerning the parameters related to the Pomeron interaction, we use the same ones given 



in Ref.@, such as F(t) = l^fi-t/artf ■ ~ m ") = ^TCTT with = 1-OGeV. The 



4m2-( (l-t/0.71) 2 ' J \f ""J fj%-(p 2 -m 2 ) 

trajectory of the Pomeron is given as 

a(t) = 1.08 + 0.251 (10) 

As for the quark-Pomeron coupling, the light quark-Pomeron coupling is fixed as (3q 2 ~ 
3.5GeV _2 |l9|. However, it is known from experiment that the effective coupling of the Pomeron 



to heavy quark is rather weaker than the one of a light quark[|24[. Therefore, we use (3 ~ 
(nT>u,d/nT>c)Po ~ 0.23/?o for the charm quark-soft Pomeron coupling, by taking account of the mass 
effect of charm quark propagator in the Pomeron-charm quark vertex. Other parameters are 
fixed as m c = 1.5 GeV and Q 2 = (2m c ) 2 . 

As mentioned above, for the polarized gluon distribution function AG(x, Q 2 ), we take two typ- 
ical parameterization models of GS96|9| and GRS V96 [|Hj] , while for the unpolarized distribution 
G(x, Q 2 ), we take the model of GRV94[Q. For the polarized and also unpolarized fragmentation 
functions, AD A + / c (£') and -D A +/ C (£'), we use the same function of the model of Peterson et al. 
taken up by the Particle Data Group in Ref . [[23] . 

Calculated cross sections and 2-spin asymmetries are presented in Figs. 2 and 3. Since 
the variables y and s are included in the amplitude mostly as a factorized form ys, the cross 
sections have almost the same behavior for the same values of ys. Since we are interested in 



the target fragmentation regions, we have shown only the case with z = 0.1, whose kinematical 
range mostly covers the region with positive rapidity for the produced A+. Note that we have 
taken the z axis to be the direction of the target proton p 2 and thus, the region with positive 
rapidity for the produced A+ corresponds to the target fragmentation region. As shown in Figs. 2 
and 3, the 2-spin asymmetry A^ C L rather largely depends on the model of AG(x, Q 2 ) in some 
kinematical regions. Therefore, the 2-spin asymmetry for the proposed process is promising for 
testing the gluon polarization. Of course, the analysis might be rather primitive because the 
present calculation is limited in the lowest order. To get more reliable predictions, in addition 
to the next-to-leading order calculation, we have to refine the Pomeron model and also to have 
a good knowledge of the polarized fragmentation function of a c quark to A+ decays. Although 
these subjects have their own interest and need further investigation, they are out of scope for 
the present work. 

In summary, we have calculated the cross section and the 2-spin asymmetry for diffractive 
A+ productions in polarized pp reactions for the target fragmentation region at y/s = 50GeV 
and 500GeV. We found that the calculated results largely depend on the model of AG(x, Q 2 ) in 
some kinematical regions and thus, the process is quite promising for extracting information on 
polarized gluons in the proton. 

Although in this work calculations were carried out expecting the forthcoming RHIC experi- 
ments, the same analysis might be applied also for the proposed HERA-N experiments. 

References 

[1] For a recent review, see B. Lampe and E. Reya, |hep-ph/9810270| (1998). 

[2] Spin Muon Collaboration, D. Adeva et al, Phys. Rev. D58, 112001 (1998). 

[3] E142 Collaboration, P. L. Anthony et al, Phys. Rev. D54, 6620 (1996). 

[4] E143 Collaboration, K. Abe et al, Phys. Rev. D58, 112003 (1998). 

[5] E154 Collaboration, K. Abe et al, Phys. Rev. Lett. 79, 26 (1997). 

[6] E155 Collaboration, P. L. Anthony et al, SLAC-PUB-8041, |iep-ex/ 9904002] (1999). 



[7] HERMES Collaboration, K. Ackerstaff et al., Phys. Lett. B404, 383 (1997); A. Airapetian 
et al, Phys. Lett. B442, 484 (1998). 

[8] R. Mertig and W. L. van Neerven, Z. Phys. C70, 637 (1996); W. Vogelsang, Phys. Rev. 
D54, 2023 (1996). 

[9] T. Gehrmann and W. J. Stirling, Phys. Rev. D53, 6100 (1996). 

[10] M. Gliick, E. Reya, M. Stratmann and W. Vogelsang, Phys. Rev. D53, 4775 (1996). 

[11] G. Altarelli, R. D. Ball, S. Forte and G. Ridolfi, Nucl. Phys. B496, 337 (1997); Acta Phys. 
Polon., B29, 1145 (1998). 

[12] C. Bourrely, F. Buccella, O. Pisanti, P. Santorelli and J. Soffer, Prog. Theor. Phys. 99, 1017 
(1998). 

[13] E. Leader, A. V. Sidorov and D. B. Stamenov, Int. J. Mod. Phys. A13, 5573 (1998); Phys. 
Rev. D58, 114028 (1998); Phys. Lett. B445, 232 (1998). 

[14] L. E. Gordon, M. Goshtasbpoul and G. P. Ramsey, Phys. Rev. D58, 094017 (1998); Asym- 
metry Analysis Collaboration, Y. Goto et al., Polarized Parton Distribution Function in the 
Nucleon, |hep-ph/0001046| (2000). 



[15] See, for example, E. Predazzi, |hep-ph/9809454l (1998). 



[16] RSC Coll., Proposal on Spin Physics using the RHIC Polarized Collider, BNL, August 1992. 



[17] W.-D. Nowak, [hep-ph /96054 1 1| , published in 'Trends in Collider Spin Physics ' Proceedings. 



by Y. Onel, N. Paver, A. Penzo. World Scientific (1997). 

[18] N. I. Kochelev, T. Morii and T. Yamanishi, Phys. Lett. B405, 168 (1997). 

[19] A. Donnachie and P. V. Landshoff, Phys. Lett. B191, 309 (1987); Nucl. Phys. B303, 634 
(1988). 

[20] C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas, Phys. Rev. D27, 105 (1983). 

[21] G. A. Schuler, Nucl. Phys. B299, 21 (1988). 

[22] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67, 433 (1995). 



[23] The Particle Data Group, Eur. Phys. Jour. C3, vol.3, 1 (1998). 
[24] A. Donnachie and P. V. Landshoff, Phys. Lett. B207, 319 (1988). 

[25] A. Donnachie and P. V. Landshoff, Nucl. Phys. B231, 189 (1984); Nucl. Phys. B267, 690 
(1986). 



p 




Figure 1: Lowest diagram for p + p — > p + A+ + X. 
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Figure 2: 2-spin asymmetry and cross sections for y/s = 50GeV and z = 0.1. Solid line is for 
GS96 and dashed for GRSV96. In the right figure, unpolarized cross section is also given in 
dotted line. 
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Figure 3: The same as in Fig. 2, but for y/s = 500GeV. 



